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Abstract 

The observed relationship between stellar mass and effective radius for early 
type galaxies, pointed out by many authors, is interpreted in the context of 
Clausius' virial maximum theory. In this view, it is strongly underlined that 
the key of the above mentioned correlation is owing to the presence of a deep 
link between cosmology and the existence of the galaxy Fundamental Plane. 
Then the ultimate meaning is: understanding visible mass - size correlation 
and/or Fundamental Plane means understanding how galaxies form. The mass 
- size relationship involves baryon (mainly stellar) mass and its typical dimen- 
sion related to the light, but it gets memory of the cosmological mass variance 
at the equivalence epoch. The reason is that the baryonic component virial- 
izes by sharing virial energy in about equal amount between baryons and dark 
matter, this sharing depending, in turn, on the steepness of the dark matter 
distribution. The general strategy consists in using the two-component tensor 
virial theorem for determining the virialized baryonic configurations. A King 
and a Zhao density profile are assumed for the inner baryonic and the outer 
dark matter component, respectively, at the end of the relaxation phase. All 
the considerations are restricted to spherical symmetry for simplicity. The effect 
of changing the dark-to-baryon mass ratio, m, is investigated inside a ACDM 
scenario. A theoretical mass - size relation is expressed for the baryonic com- 
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ponent, which fits fairly well to the data from a recently studied galaxy sample. 
Finally, the play of intrinsic dispersion on the mass ratio, m, is discussed in the 
light of the cusp/core problem and some consequences are speculated about the 
existence of a limit to/ expected by the theory. 
Keywords: 

Celestial Mechanics, Stellar Dynamics; Galaxies: Clusters. 



1. Introduction 

The physical grounds on the basis of Clausius' virial maximum theory (CVMT) 
lies in that: to explicitely consider the additional degree of freedom the problem 
has, owing to the existence of one further component really present in a galaxy 
structure: the dark matter halo. This is not to be included into the system in 
a way to handle it as a single one and then applying the one-component tensor 
virial theorem. On the contrary, the system is splitted into two subsystems: 
one of baryons, the other one of non-baryonic dark matter (DM) to which the 
two-component tensor virial theorem is applied. Then the double system gains 
some new ways to share out the amount of potential energy which has to ap- 
pear in the two virial equilibrium equations. Briefly summarizing, the general 
strategy consists to use the two-component tensor virial theorem (e.g., Brosche 
et al., 1983; Caimmi & Secco, 1992) to describe the virial configuration of the 
baryonic component embedded in a DM halo at the end of relaxation phase 
(see, Bindoni & Secco, 2008). It reads: 

2(T u )y = (K)«; {u = B,D;i,j = x,y,z) (1) 

According to the scalar virial for one component, the potential energy tensor, 
which has to enter into the tensor virial equations, is the Clausius' virial (CV) 
tensor, (V u )ij, built-up of the self potential-energy tensor, (fi u )jj, and the tidal 
potential- energy tensor, (V uv )ij (u = B , D; v = D,B). Then, according to the 
scalar virial theorem, the trace of CV tensor, related to the inner bright, B, 
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component, has to be read: 



v B = n B + V BD 



(2) 




- dxi = J p B {r~B ■ /s) dxi 
dxB = / p B (r B ■ Id) dxh; 



(3) 



(4) 



where ps is the B component density and f B , Jd are the force per unit mass 
due to the self and DM gravity, respectively, at the point r B and $b, &d 
are the related potentials. Conversely, the total potential energy tensor of the 
B component is: (Q B )ij + (W B D)ij, where the interaction energy tensor is: 

ij dx B ; and the potential tensor (e.g., Chandrasekhar, 



To be noted that in general: (V B o)ij ^ {W B o)ij, the difference gives the 
residual energy tensor (Caimmi & Secco, 1992). 

In spite of the total potential energy: (E pot ) B — Q B + W B d, which has an 
increasing monotonic trend by increasing the ratio between the volume taken 
by baryons over that of DM halo (Secco, 2005, hereafter quoted as LS5, Fig. 
2), the trace of CV tensor, V B , shows, under some constraints, a maximum 
(CVM) for a special virial configuration. That means, in turn, a minimum of 
the random kinetic energy to obtain equilibrium for the baryonic component 
when it is completely embedded into a DM halo assumed with a fixed density 
profile and dimension. The possibility to get a maximum is subjected to two 
restrictions, namely (1) the DM density distribution is decreasing not too fast 
inside the bulk of baryonic matter, and (2) the total DM mass maintains above 
a threshold. Indeed we are dealing with an amount of CV energy to be divided 
between two components: that in baryons and that in DM fraction inside the 
baryon container (due to Newton's first theorem). If a power- law DM density, 
Pd ~ l/ r ' increases too fast toward the center, the weight of tidal energy 
tends to overcome that of the self gravitational energy due to baryons, until a 
divergence in the energy ratio occurs as soon as d — > 2 (Marmo & Secco, 2003). 



1969) due to the DM is: ($ D )ij = GJ p D {x') 



(x 




D ■ 
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To explain some of the main Fundamental Plane (FP) features in the light of 
the CVMT theory, a DM distribution characterized by d ~ 0.5 is needed (e.g., 
Bindoni, 2008). 

Moreover, decreasing the DM total amount the sharing of virial energy needs 
that the DM fraction inside baryon container has to become all the available DM 
amount. In other words the radius corresponding to the CVM moves outwards 
so that the two volumes coincide. Under this thresold CVM doesn't exist. 

In sect. 2 the existence of a special tidal dimension for the baryonic compo- 
nent due to the gravitational interaction of DM is underlined and a model to 
link it to the light dimension r e is described. At sect. 3 the ACDM cosmological 
scenario is introduced and the typical CVMT scaling laws of r e , I e , <r with the 
only two free parameters Mb, m are considered by means of the only two expo- 
nents: the slope of DM density profile and the cosmological local slope of mass 
variance. At sect. 4 and 5 theoretical fits vs. observations are taken into ac- 
count with some linearization procedure and some useful relationships. In sect. 
6 the theoretical slopes of mass-size relationship are successfully compared with 
those Tortora et al. (2009, hereafter quoted as Tal09) obtained by observations. 
The role of the intrinsic dispersion due to m as function of the slope of DM 
density profile are discussed and compared with CVMT's expectations for some 
of the main scaling laws, in sects. 7,8,9. The discussion and conclusions follow 
at Section 10 and 11, respectively. 



2. Sharing virial energy between interacting components 

It is not surprising that a relationship between the visible mass, mJ* and 
its spatial dimension does exist for early type galaxies. Indeed, according to the 
CVM theory's expectation a special scale lenght is induced, when the maximum 
exists (see, Secco & Bindoni, 2009, hereafter quoted as SB9), from the dark 
matter halo on the baryonic gravitational field. In this way the gravity, which 



1 Mostly in stars, assumed to form very early. So the visible mass Mb — M* 
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for its nature has an infinity range, acquires an intrinsic macroscopic scale. 
Among the infinite possible dimensions corresponding to virial equilibrium of a 
B component embedded inside a dark one, D, only one is able to maximize the 
Clausius' virial energy by sharing it in about equal amount between baryons 
and dark matter particles. The consequence is that visible virialized matter 
does not have a whatever dimension but a special one, that is the tidal radius 
at, strictly related to that of DM virial radius, ajj. 

The assumed DM density profile is a special case of Zhao (1996) family (for 
more insight see, Caimmi et al., 2005): 

where r a u is a scale radius and p o — p{t d)- More specifically, the profile of 
Eq. ([5]) has a central core and tends to a power-law density profile with exponent 
d for sufficiently large radial distances. We will come back to the relevance of 
the d-value in the next sub-section. 

The related DM virial radius - tidal B-radius relation is found to be (SB9): 

(vnB 1 Mb \ 3=2 
at = M~ D J a ° (6) 

where a t , ap, are the tidal and the virial radius, respectively, M is the total 
mass, uqb the self-mass distribution coefficient, v' v the reduced interaction coef- 
ficient, and the indices, B, D, denote the star and DM subsystem, respectively. 
The reduced interaction coefficient is defined as: 

, , v v v (x) a B M D 

rai j a ao Mb 

where vy is the interaction coefficient, x the B to D virial radius ratio, and m 
the D to B total mass ratio. It can be seen that v' v (x) w const to a first extent 
(SB9). For deeper insight to Eq. ([5]), further considerations are needed. A 
dimensionless trace of CV tensor can be obtained dividing both sides of Eq. ([2]) 
by a normalization energy, GM^F/ao, where G is the gravitational constant 
and F a form factor (F=2 in the case under discussion of spherical-symmetric 
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configurations). The result is: 

y B= _^_VjM (8) 

X X 

The condition: — o, yields the value at expressed by Eq. ([6]). It is note- 
worthy that existence of at which means from one side the existence of the 
maximum in CV energy and then its sharing between the two components, im- 
plies from the other side the existence of a minimum value for to. Indeed, if the 
baryonic component has to be embedded inside the DM halo (according to the 
tensor virial theorem extended to two componens) the location limit of CVM 
is the DM component border. Studing how changes the ratio Xt = at/ an as to 
changes it is easy to prove that Xt grows at decreasing to so that at reaches ao 
for the lowest value: 

corresponding to Xt = 1 in Eq. ([5]). That in turn leads to the expectation of 
zones of exclusion in the scaling laws based on the CVM presence. In fact, 
values to < to; would imply a DM halo embedded within the star subsystem, 
contrary to both current cosmological scenarios and observations such as flat 
rotation curves in galactic disks which extend well outside the visible region. 

It is to underline the strong analogy which exists between the tidal radius at 
and that von Hoerner (1958) found for a spherical star cluster embedded in the 
Galaxy tidal potential. At first sight the result might appear surprising because 
in the cluster case we are dealing with two forces, the self gravitational- and 
tidal-force which are acting in opposite directions; here, on the contrary, the 
tidal radius is a consequence of two attractive forces toward the same central 
direction: the self gravity per unit mass fs and the other one /d due to DM 
both applied to each point tb fEqs l3l4l) . But the balance which defines at is 
not between forces but between energies, given by forces times positions: self 
and tidal. They have different trends as as decreases. 

Tidal energy of B due to subsystem D, goes as follows: according to first 
Newtons' theorem the fraction of D mass which exerts dynamical effect on the 
B component is only that inside the B boundary, which for spherical-symmetric 
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configurations is proportional to Mi>(^B-) 3 ~ d , then the integral of force times the 
position decreases toward like Mr>(^) 2 at decreasing a^. The other limit 
that is the maximum value of such an energy, corresponds to the configuration 
in which the whole D mass is inside the B component i.e. the two boundaries 
coincide. Exactly the contrary occurs to the self energy: it grows as as decreases 
and tends to zero at increasing as- 

Then in a configuration where as > at, the tidal energy term in the 
Clausius' virial (EqJ2]) is dominant in respect to the self-energy one. When as 
becomes less than a t , it occurs the overwhelming of the self-energy over tidal- 
energy. All that turns to be in strict analogy not only with von Hoerner's radius 
but also with the Hill's radius (see, Binney & Tremaine, Chapt.8, 2008). Indeed 
on the basis of both definitions it lies the restricted three body problem. But it 
is to be noted that our new definition of tidal radius (Eq|6]) refers not to two off 
centre objects but to two concentric ones. Then it looks like a generalization of 
those previously given. 

2.1. CVMT and DM cusp/core problem 

Whithout entering the complex, still open problem (see, e.g., Bindoni, 2008) 
we underline only the relevance it has inside the CVMT. With regard to DM 
density profiles characterized by the presence of a central cusp, such as NFW 
(Navarro et al., 1996), the term "cored" has to be intended as related to milder 
(pDM r~ x l d , d < 1) slopes in the central region instead of erasing the cusp. 
A similar situation occurs in the case under consideration. With regard to 
DM density profiles characterized by the presence of a central core, such as 
the one expressed by Eq.©, a "cuspy" behaviour is shown by steeper (d > 1) 
slopes at sufficiently large (r >> r £>) radial distances without erasing the core. 
What is relevant is indeed the slope of DM inside which the bulk of baryonic 
matter in stars is embedded because the gravitational interaction may lead for 
the two cases: d < 1, d > 1, to completely different results. Coming back to 
the normalized virial energy (j8|) and getting advantage by the approximation 
v' v const, see Eq. ((TJ), the ratio £t between the self and the tidal energy at x t , 
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becomes (Marmo & Secco, 2003): 

C t ="*~(2-d) (10) 
Vbd 

which means the exact balance between self and tidal contibution is reached 
when the logarithmic slope of the DM density profile is d = 1 (i.e., cuspy 
profile). Then Vbd/^b increases towards oo as soon as d — > 2. Moreover, 
it is easy to prove that in this cuspy case, d=l, the contribution of DM tidal 
component to virial energy becomes of the same amount of the baryonic one, 
whichever is the ratio m. In other words this case is equivalent to handle with a 
unique component in the virial equilibrium equation having an equivalent self- 
mass distribution coefficient, Dq = 2 • z^ns, two times that of the single baryonic 
one. It means that a double system considered on the maximum of its inner 
component virial energy may degenerate into a single system only in this sense: 
not if DM mass reduces to (because the DM amount may be reduced only 
until to mi) but as soon as the DM component is able to share exactly with the 
inner one the amount of its virial energy. 

2.2. The r e — at link 

The relationship r e — M* can be obtained via Eq.©, which links the size ratio 
between the two components to their mass ratio, provided a connexion exists 
between r e and the tidal radius at, by a model, within a cosmological scenario 
where the DM halo dimension is related to the density primordial perturbation 
spectrum. 

We assume that the ETGs are members of a perfect homologous family and 
then for each galaxy we handle with three different profiles for baryonic, light 
and DM distributions, each of them assumed to be the same for the whole family. 
The universality of our choice might be questionable in particular that of DM 
either for the still pending cusp/core problem or for the difficulty to reformulate 
the CV theory when baryons are embedded inside a halo density profile very 
different from that of a power-law. 
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Here we refer to our non-linear spherical model considered in SB9 where it 
was assumed for the bright B-component the empirical surface light density law 
proposed by King (1962) as: 

I{R) = kh { [l + {R/Rc) 2 ] 1 ' 2 ~ [1 + {Rt/Rc) 2 ] 1 ' 2 } (U) 
where Rt is the value of projected radiuj^ R at which / reaches zero. The law 
has the advantage to take into account the existence of a cut-off in the surface 
light distribution as one expects for globular clusters (GC) (the profile (TlTj) was 
born indeed for GCs and then it has been extended to ETGs). As shown in 
previous subsections, that is particular suitable because we expect for ETGs the 
presence of a cut-off which is at and has an analogous role of von Hoerner's one. 

The value R c corresponds to the core-radius and the fcz, value is linked to 
the central surface light density I a by: 

'■^'('-Fskfpf (12) 

The corresponding baryonic matter projected density is: 

m) = kM " [i + (im 2 ) (13) 

which is linked to the spatial baryonic matter density by the Abel integral 
equation (Binney & Tremaine, 1987, Chapt.4): 

It should be noted that in our model kM 7^ being one of our main as- 
sumptions (see, SB9). Integration of S(i?) with respect of 2irRdR gives the 
total projected mass within the projected distance R from the center which 
becomes the luminosity function, L(X), with the substitution of S(i?) — > I(R) 
and kM — > fc^: 

L(X)=Trr 2 c k L F L (X) (15) 



2 Due to the choice of spherical model the projected radial quantities: Rt, R c , Re, coincide 
with the corresponding three dimensional ones: r t ,r c ,r e . 
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where: 

(l+X) 1/2 -l X 



Fl{X) 



ln(l + X) -4- 



{l + X t ) l/2 1 + Xt 



(16) 



-o 2 • *>=m ™ 

According to Eq. (fT5| . as At >> 1 the limit of L{Xt) goes approximately to: 

L{X t )^nrlk L ( 18 ) 

The square of the effective radius r e normalized to r c is given as solution of the 
equation: 

L(X e ) = \r 2 c k L F L {X t ) (19) 

which links also r e to the cut-off radius Rt — rt — at- 

Referring to DM distribution its power-law is given by Eq.([S]) so that in 
normalized form it becomes of this kind: 

/d(£d) " TTTcb? (20) 

where, £d = t/o>Di C'd = o-bItoD is the DM concentration, d < 1 and the 
normalization is done to the density value p {£,D = 1/Cn) at the scale radius 
r D- As soon as the concentration of the King's component is fixed for the whole 
galaxy famihjfl, then at turns to be linearly proportional to r e . Now the next 
step regards the key point that is to connect an of Eq. © to the DM halo mass 
by cosmology. 

3. Cosmological framework 

It is worth recalling that the key point to explain the Mass — Radius rela- 
tionship comes out from the cosmology. If otherwise, neither all the well known 
scaling relationships for galaxies nor the FP, which implies their existence, could 



3 According to SB9, if wc assume: rt/r c = 10 that is the B-component concentration, 
Cb = 10, then it turns to be: r e /r c = 1.70. 
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be interpreted (see, LS5, SB9). That was one of our first and most important 
consequences of CVMT, in spite of the apparent paradox that the FP, as a 
whole, does contain a degeneracy in respect to the initial density perturbation 
spectrum (first pointed out by Djorgovski, 1992). 

For a scale-free power spectrum, P(k) <~ K n , the mass variance of DM be- 
comes: 

4(t)~fl 2 (t)M-(" +3 »/ 3 (21) 

where the self-similar linear growing factor, D(t), allow us to describe the gen- 
eral density perturbation as: 5(x, t) = 5(x)D{t). 

3.1. ACDM scenario 

We will refer to the cosmological standard model derived from WMAP pre- 
cision data only (Spergel et al., 2003; Spergel et al., 2007; Binney & Tremaine, 
Chapt.9, 2008) for a flat A-dominated universe: i.e., ACDM model defined by 
the following parameters: <r% = 0.9 ± 0.1; h = 0.72 ± 0.05; matter density, 
n m h 2 = 0.14 ± 0.02, baryon density, Vt b h 2 = 0.024 ± 0.001 and primordial 
spectal index n s = 0.99 ± 0.04 which is consistent with the Harrison-Zel'dovich 
scale- invariant value (n s = 1). The trend of mass- variance spectrum o~m d at 
equivalence between matter and radiation (after that epoch the microphysics 
may never affect the DM fluctuations) is shown in Fig. 1 as function of Mo- 

In this self- similar model, an effective final spectrum index, n e ff = n e , may 
be obtained by local slope of ctm d at a given Mu value, as follows: 

- = -^I^ = <". + ^ <»> 

At the halo virialization the following scaling laws hold (see, SB9): 

a D ~M D h '; (23) 

Pod - M D ^ +3)/2 (24) 
1/ y = (5 + "e) = 3a eg + l (25) 

where an is the halo virial radius and p Q u = p(r D) according to an assumed 
Zhao's (1996) density profile. 
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Figure 1 : The trend of mass variance <tm d i at equivalence epoch, as function of Mo , for a 
flat A— dominated universe derived from WMAP precision data (ACDM model defined by the 
following parameters: ct 8 = 0.9±0.1;fc = 0.72±0.05; Q m h 2 = 0.14±0.02; fl b h 2 = 0.024±0.001; 
n a = 0.99 ± 0.04; see text). 
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Figure 2: The trend of the local slope, a eq , Eq. (1221 (top panel) and of 7', Eq. I|25jl (bottom 
panel) as function of DM halo mass in the mass range here considered. 
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3.2. Scaling laws 

Coming back to Eq. © the scaling law of r e with Mb for the whole homology 
family reads: 



a * ~ \Tr) a ° , — ~ m '' M b ( 26 ) 

owing to (|23[) . by definition: Mjj — wiMb and under the assumption that bright 
and dark matter profiles are universal {y' v and vqb are depending on both and 
inner profiles, respectively). It is to be underlined that in the CVMT all the 
three main quantities of FP: r e , I e , <j , may be expressed in terms of only two 
independent parameters, as Brosche as pioneer discovered (Brosche, 1973), e.g. 
m and Mb (LS5, SB9), and have to scale with them by means of exponents 
in which enter only two quantities: d, i.e., the densiy profile of DM and 7' 
related to the cosmological local slope of mass variance. Indeed the following 
nice scaling laws have to hold (LS5, SB9): 



3 - d 



m l M B ; I = i = 2 7 ' (3 f = -2r (28) 
7 (3 - d) 



s 1 (3 - d) - j iy-i 

a ~mM B ;s=-~ — _ rf) ; S = - — 
To be noted that the bright density distribution doesn't enter. 



(29) 



4. Theoretical fits to observations 

The best pioneeristic attempts to find a relationship M* — r e appear in the 
literature due to Chiosi et al. (1998) and Tantalo et al. (1998) starting from 
the data of Carollo et al. (1993) (see, Contini, 2008). It reads as follows: 

r e = 17.13 (Af B , 12 ) a557 (30) 

where r e is in kpc and M B ,i2 = M B /(10 12 M ) . The data refer to 42 ETGs, 11 
SO + 1 which probably could be a spiral. The fit is given in Tantalo et al. (1998). 
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As we will see, relationships of kind (|5P|) may gain a full explanation inside the 
framework of CVMT. We refer here to effective radii and star subsystem masses 
determined in a recent attempt (Tal09) for local elliptical galaxies, using model 
populations and a large uniform data set from a well populated homogeneous 
sample ( Prugniel & Simien (1996 (PS96); 1997 (PS97)). More specifically, 
Tal09 aim to carefully consider all the factors which may enter the FP. Many 
accurate, interesting fits and mass models are given therein. 

Better and more extended samples of ETGs are available now in the liter- 
ature, in particular the data from SDSS, where a similar mass-size relation is 
shown (e.g., Bernardi et al. 2003) while mass determination is unnecessary in 
testing the Faber- Jackson relation (e.g., Nigoche-Netro et al. 2010). For this 
reason, as a first step Tal09 results shall be used in the following. 

The related mass-size relation (Tal09, Fig. 4) when binned averaged points 
are considered, has a best fitting mean slope over the full range of the data, equal 
to: 0.58±0.05, in agreement with previous estimates of typically ~ 0.6 (Bernardi 
et al. 2003; Shen et al. 2003; Mamon & Locas 2005; Napolitano et al. 2005) and 
recent estimates of typically ~ 0.55 (Bernardi et al., 2011, hereafter Bernll). 
Moreover two mass regimes are identified: faint galaxies (M* < 10 11 ' 1 Mq) with 
a slope: 0.36 ± 0.13, and bright galaxies (M* > 10 111 Mq) with a greater slope: 
0.73 ±0.12. At fixed mass ratio to, this trend appears immediately in agreement 
with the theoretical trend of I/7' expected by the scaling relation (f2~T|) where the 
exponent increases at increasing Mo because the local slope of the cosmological 
mass variance increases (Eqs. (|22|) . ([25"])), 

More precisely the CVMT is able to differentiate the contribution of either 
mass ratio to or Mb (see, (|27j) ). both having a variable exponents in order to 
produce a continuous set of curves on the data plane (log M*,log r e ) (Tal09, 
Fig. 4). That allows us to improve the resolution of the relationship we are 
looking for, taking into account the trends of the points which show, in Tal09's 
data, the individual galaxies of the sample. According to the theoretical scaling 
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Figure 3: Mass-Radius data of local ETGs from PS96 (see, Tal09). The solid lines in colours 
are the theoretical curves given by Eg . (13 IB when the mass ratio m assumes (bottom- up) the 
values: 1,4,10, and d = 0.5. The red square refers to the calibration point for which it 
is assumed m° = 4 so that the other corresponding parameters entering 1311 1. are: (M B = 
1 ■ lO n M , r° = 3.133 kpc, r a = 0.076, Y Q = 2.103). The bended dashed line is the fit given 
by Tal09, corresponding to the two regimes of low and high masses. The dot-dashed straight 
line shows the mean mass-size relation found by the same authors. On the top left corner the 
rms errors are shown as a cross. The limit case m = 0.5 ~ m; is also given as lowest curve 
delimiting the underlying avoidance region. See the on-line version for colours. 
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- = Km r M^ (31) 



law (|27|) . we have: 

_ m r Mj 
( m ) (Afg)'- 

if = r° -j— = const. (32) 

(m°) ro (M»)^ 

The values: (Mg, r°, r Q , j' ) refer to an arbitrary calibration galaxy chosen to 
normalize the trends in the sample. 

4-1. Calibration 

Because in Tal09 M B ~ 10 11 M Q works as characteristic mass scale between 
the two mass regimes of faint and bright subsamples, we select it as mass cal- 
ibration point. The corresponding effective radius is r° = 3.133fcpc which has 
to read on the mean mass-size relation found by Tal09 with the mean slope 
ctM = 0.58 ± 0.05. About the mass ratio m°, there is not an a priori value to 
be chosen then it will be considered a free parameter even if some hints may 
come from cosmology. It is to be underlined that the values of both r 0l ^' are 
determined by m° because it settles the amount of dark matter halo which in 
turn gives the local slope ([22)1 in the cosmological mass variance ctm d ■ 

According to the mean 68% confidence range given by Spergel et al. (2003) 
for matter density and baryon density in their ACDM model, we obtain a cor- 
responding range for the cosmological ratio: Md/Mb as: 

Q b h 2 



1 



m = "™ fe2 ; 3.8<m<5.9 (33) 

Q m h 2 

As scouting sake a more extended range (m = l-=-10) will be considered assuming 
as reference value: m = 4. Moreover, according with the CVMT we know that 
an inferior m— limit has to exist, given by Eq.(|9"j), which for the cases here 
considered turns to be in the range: 1 ~ 0.5 (Bindoni, 2008). For sake of 
investigation the range assumed will be: m = 1 -j- 10. For special aims we will 
consider also the extreme limit m = 0.5 ~ mi. 
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Jj.,2. Theoretical relationships 

Three mass-size theoretical relationships of kind (|31l) may be performed by 
assuming different values of m° for the calibration point in each of them (e.g., 
m°=j=l, 4, 10). Starting from each of these curves with label j Eg. ([31]) yields 
a continuum of curves by changing m. For sake of simplicity we consider only 
a couple of them for each j: m=4, 10; m=l, 10; m=l, 4, respectively. In Fig. 
3 the triplet of curves corresponding to: m°=j=4; m=l, 10 is shown. They are 
drawn at fixed m but with the exponents r and I/7' changing with Mb- The 
constants (f3"2|) are: 



The common property of all the curves is to have a positive local slope increasing 
with m at fixed value of Mb- The triplet of curves, shifts towards bottom as 
soon as the value of m° attributed to the calibration point increases from 1 to 
10. The case m° = 4 seems to produce a better fit to the observed spread which 
remains greater than the typical vertical rms error bar shown in Fig-El 

4-3. Linearization 

Our aim is to linearize the theoretical curves (|35[) in order to obtain useful 
numerical relationships. In logarithmic form the set of equations becomes: 



K* = K3(r° el m° = 3, M° B ,r , j' ); j = 1,4, 10 



(34) 



The equations of the corresponding curves become (i=l, 4, 10): 



K 1 (mlY Mj ;m° = 3 = 1; m\ = 1, 4, 10 
< K 4 (mi) 1 " Mj ;to° = j = 4; m\ = 1, 4, 10 
K io ( m jy M J ]m ° =j = 10; mf = 1, 4, 10 



(35) 



V, 



i (x) = Ki + H (x) ml + (j (x) x- (*, j = 1, 4, 10) 



(36) 
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with the following definitions: 



(37) 



V = log r e 
x = log Mb 
${x)=r{M B ) 

C'W = 7' (M B ) 

ml — log m\ 
Ri = log Ki 

They may be linearized by fixing a reference point in each curve with coordinates 
B i^i > Vi ) an d taking the following local slope at the same point : 



pUx) 




x + ciix) 



(38) 



It is noteworthy that the theoretical slope at a reference point x preserves mem- 
ory of the cosmological mass variance either through its first derivative (see, 
Eqs. (|22| [25 ]) ). inside the third contribution (£ term in Eq. (|38"1) ) or its second 
derivative (first and second term of the same equation) . The transformed linear 
equations of set ([33)1 become: 



z\{x) = K' i +pl(x) x 



- V'3 



K J M 



(39) 



where: 



K' J - - 



3 ^3 



z{ (x) 



(40) 
(41) 



It may be interesting to distinguish between the two contributions: one due to 
the derivative of 1/j'(Mb) and the other one to the derivative of t(Mb)- Both 



In our curves the x-coordinate always corresponds to x 3 . = x = 11.00. 
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added to the term £| (x) enter the slopes p\(x). To this aim we define: 




^(5) = ^(x)+^'(5) + Cf(5) (44) 

The different contributions are evaluated in Tab. 4. It is noteworthy that the 
trend of £(x) = y^p^ji which determines both the slopes q(x), g(x) and in 
turn p(x), is depending on the value of Mjj. That means: at fixed Mb, the 
slopes have to be the same if to is the same. So curves with the same label i 
and different j are simply the same curves vertically translated one in respect 
to the other owing to the different values of log r e . Furthermore, at fixed j, 
a curve at higher to is simply obliqually shifted over that at lower to, towards 
higher values of Mb- Moreover, at fixed d, whichever it is, by definition of (l27|) , 
it occurs: 

That allows us to collect in a unique Tab.l the trends of 7', dt^/dx and p in 
the three cases to = 1, 4, 10 instead of handling with three tables. 

5. Useful Relationships 

According to the linearization procedure performed in subsect.4.3, we may 
relax the condition of variability of r and 7' along the Mb coordinate, as occurs 
in the theoretical Eqs. (|35| . and consider the three straight lines passing through 
the calibration point having, at this point, the same slopes of the three curves 
of the set (j35|) with: j = l,m = 1; j = 4, m = 4,j = 10, to = 10. The linear 
equations, as approximations of Eqs. (|35j) . are, respectively: 



kpc 



15.240 (M Si i 2 ) - 687 ; to = 1 

17.824 (M Bil2 ) - 755 ; to = 4 (46) 
20.137 (M B , i 2 ) - 808 ; to = 10 
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Table 1: The trends of 7' and the first derivative of ( = I/7' as function of the baryonic mass 
Mg and m are shown. The end product p, according to 1441 . is also given for any d (see, 
Figs.3,5). 



log M B 


Y(x) 


7 '(x) 


i{x) 


dQ{x)/dx 


d((x)/dx 


dQ(x)/dx 


p(x) 


p(x) 


p(x) 


X 


m=l 


m=4 


m=10 


m=l 


m=4 


m=10 


iii=l 


m=4 


m=10 


9.9 


2.262 


2.210 


2.174 














10.0 








0.016 


0.018 


0.020 


0.603 


0.650 


0.687 


10.1 


2.245 


2.192 


2.155 














10.2 








0.017 


0.019 


0.022 


0.623 


0.668 


0.715 


10.3 


2.228 


2.174 


2.134 














10.4 








0.018 


0.021 


0.023 


0.636 


0.694 


0.735 


10.5 


2.211 


2.154 


2.113 














10.6 








0.019 


0.022 


0.025 


0.653 


0.712 


0.763 


10.7 


2.193 


2.134 


2.091 














10.8 








0.020 


0.023 


0.025 


0.671 


0.739 


0.781 


10.9 


2.174 


2.113 


2.069 














11.0 








0.020 


0.024 


0.027 


0.687 


0.755 


0.808 


11.1 


2.155 


2.091 


2.047 














11.2 








0.022 


0.025 


0.029 


0.715 


0.781 


0.843 


11.3 


2.134 


2.069 


2.023 














11.4 








0.023 


0.027 


0.031 


0.735 


0.813 


0.879 


11.5 


2.113 


2.046 


1.998 














11.6 








0.025 


0.029 


0.033 


0.763 


0.847 


0.915 


11.7 


2.091 


2.022 


1.972 














11.8 








0.025 


0.031 


0.034 


0.781 


0.880 


0.941 


11.9 


2.069 


1.997 


1.946 














12.0 








0.027 


0.032 


0.036 


0.808 


0.905 


0.980 


12.1 


2.047 


1.972 


1.920 














12.2 








0.029 


0.034 


0.038 


0.843 


0.941 


1.028 


12.3 


2.023 


1.946 


1.892 
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Table 2: Relative scatter between values from the straight lines of Eqs. (14611 and 
theoretical curves of Eqs. 1351 . at fixed coordinate log Mg and mass ratio m. 



from the 



log Mb 


m 


r* (kpc) 


r e (kpc) 




9.90 


1 


0.617 


0.550 


11 


12.30 


1 


31.189 


24.493 


21 


9.90 


4 


0.531 


0.463 


13 


12.30 


4 


40.137 


30.023 


25 


9.90 


10 


0.476 


0.405 


15 


12.30 


10 


49.333 


35.184 


29 



where Mb, 12 = Mb/(10 12 Mq) and d = 0.5. It is to be noted that the second 
equation is not too far from the preliminary pioneeristic one found by Tantalo et 
al. (1998) using the data of Carollo et al. (1993). In Tab. 2 the scatters between 
these linear approximations and the right theoretical curves are summarized. 
Then the method may be generalized thanks to the remarkable note that the 
slope, at fixed Mb, is depending only by m. So the following three general 
linear relationships may be built up: 

[y - pi (x) (x)]+pi(x) log Mb ; m = 1 
log r e = I [y - P4 (x)(x)} +p 4 (x) log Mb] m = 4 (47) 
k [y - Pw(x)x)] + pio(x) log Mb', m = 10 

which hold whatever is the calibration point (x,y), reading the corresponding 
slope, at fixed m, in Tab.l. As soon as the calibration point is (x = 11.00, y — 
0.496) the set flSJ) transforms into that of Eqs. (|46|) . 



6. Theoretical slopes vs. observations 

Tab. 3 allows us to compare our theoretical mean slopes as function of Mb 
at different mass ratio m, with those given by Tal09 in the two mass regimes, 
separated by Mb = 10 11 Mq, which are labelled with the indices: I and h 
respectively (columns (2), (3), (4) to compare with column (8)). Also the the- 
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Table 3: Mean theoretical slope values as function of Mb (column (1)) at lower (Mb < 10 11 
Mq) and higher (Mb > 10 Mq) mass regimes (index I and h, respectively), in the cases: 
m = 1,4, 10 (columns: (2), (3), (4)) are shown. The respective mean theoretical slopes in 
the whole range are given in columns: (5), (6), (7). Columns (8), (9), refer those of Tal09's 
for comparison. 



log Mb 
(1) 


(pi)i 


(pi)i 




Pi 

(5) 


Pi 

(6) 


Pw 
(7) 


T 
Pi 


f 
(9) 


(2) 


(pi)h 
(3) 


(pio)h 

(4) 


vj 
(8) 


9.90 


0.637 


0.693 


0.736 








0.36 ±0.13 




11.00 


0.762 


0.846 


0.913 


0.700 


0.769 


0.825 


0.73 ±0.12 


0.58 ±0.05 


12.30 



















Table 4: Values related to calibration red point: M^ = 1.00 ■ 10 11 Mq, r°=3.133 kpc, entering 
1131 1 and J35 D ■ The corresponding values of constants K J in i l'ill i and (135 1 > with Mb in Mq 
and K'J , which enter the numerical relationships J39H and i ]4Gll , are also given in the case 
d=0.5 when M B is in units of 10 12 Mq. 



j = m° 


j 

mi 


To 


io 


IK 


iP'(-10 5 ) 




pl(x) 


4( x ) 




1.0 


1 


0.062 


2.164 


0.462 


2.594 


15.240 


0.687 


0.0000 


0.225 


V 


4 




11 


11 


11 










J5 


10 




")■) 


11 


1-) 










4.0 


1 


0.076 


2.103 


0.476 


1.637 










V 


4 




11 


11 


11 


17.824 


0.755 


0.0145 


0.264 


J5 


10 




11 


11 


1-) 










10.0 


1 


0.086 


2.058 


0.486 


1.159 












4 




11 


11 


11 










5! 


10 


•>■> 


11 


11 


1-) 


20.137 


0.808 


0.0270 


0.297 
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oretical mean slopes over the total range are given in columns: (5), (6), (7), to 
be compared with column (9) (Tal09's data). It appears that the mean the- 
oretical slope at high mass regime is in fairly well agreement with the value 
found by Tal09 if we exclude the case of m = 10 which leads to a value out of 
the error bar. According to the ACDM scenario, the typical value expected for 
the Dark/Bright mass ratio in galaxy sample considered doesn't be around 10 
but rather about: m ~ 4 — 5. Then the theoretical case m=10 is indeed to be 
ruled out. At lower regimes the mean slope of CVMT is too high. A possible 
explanation may be the following: the mean dashed straight line of Tal09 goes 
across the theoretical lines characterized by a hxed m value (see, Fig. 3) so that 
the fit results flatter. The consequence is that also the mean slopes in the whole 
range: Pi,P4, of Tab. 3, turns to be out of the maximum limit obtained from 
observations (column (9) and dot-dashed straight line plotted in Fig. [3]). 

7. Intrinsic Dispersion vs. DM distribution 

According to the theoretical interpretation given by TCV, the (M* — r e ) 
relationship is depending on a third parameter, m, which enters pip with a 
non-irrelevant exponent r. By its definition (|27j) it is manifest that it is a 
function of n e via 7' and of d which determines the DM distribution. For a 
typical mass of DM halo on the galaxy scale, the value of n e ~ —2, then it 
follows that 7' ~ 2 and then r decreases from r ~ 1.7 at d = 0, towards r = 
when d 1, reaching the value r = 0.1 at d — 0.5. In turn it means that 
the two exponents of the (M» — r e ) relation are comparable for d < 0.5 and 
then both are to be considered. The effect of m is to shift the curves by a 
non-trivial quantity due to the role of r. To a value of .M* in each figure, it 
has to correspond an intrinsic spread in r e due to the different possible values 
of m apart from the measurement rms errors (here assumed to be of ±20%, i.e. 
±a(log r e ) = ±cr e = 0.1). It means that the same baryonic mass may have DM 
haloes of different mass values M& . 

If we move from a fixed curve of Eq. (|3~Tj) (e.g., m° — 4) this intrinsic dis- 
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Table 5: Values of the intrinsic dispersion ]49H wheighted by the two terms in d and 7^ rig- 
orously calculated as function of d starting from the case : m° = 4, 7' = 2.103 ( r = for 
d = 0.897 instead of d = 1 when 7' = 2). The up and down values are obtained for variation 
of m from m° = 4 up to 10 and down to 1. 



d 


1 

3-d 


1 

I'n 


(Alogr e ) u P 


(Alogr e ) down 


(Alogr e ) mean 


0.000 


-0.334 


0.476 


0.06 


0.13 


0.10 


0.500 


-0.400 


0.476 


0.03 


0.07 


0.05 


0.897 


-0.476 


0.476 


0.00 


0.00 


0.00 



persion, at fixed baryonic mass (e.g. M B = lO n M0) and at fixed value of r a 
{d = 0.5; i Q = 2.103, Tab.4), is given by: 

f Sr e \ Sm (3 - d Q ) - j Q 

kL= r °^ ;ro= V o (3-40 = - 076 (48) 



This scatter, for the same 7^ and the same variation of m depends on the 
exponent d tuning the DM distribution inside which the bulk of baryonic matter 
in stars is embedded. Eq.(j48|) reads as: 

( Sr e \ 1 1 5m 

In Tab. 5 the whole wheighted contribution of f4"9")l as Alogr e is shown in the 
cases with different d and changing m from m° up to 10 and down to 1 together 
with the mean dispersion scatter. It should be noted that the intrinsic disper- 
sion becomes comparable with the rms errors as soon as < d < 0.5. The 
Eq.(|49p is a differential relationship obtained by deriving the second equation 
of set ([35"1) at the mass and at fixed 7'. To underline the effect due to d, the 
same procedure may be adopted for each Mb along the curves: m — 10, 1, 0.5 of 
Fig. 3 by freezing the 7' distribution to that of the reference curve m = m° = 4 
and changing only d from to the value for which at M B (i.e., j = 2.103) r 
becomes equal (d — 0.897). The results are shown in Fig. 4. However, by 
relaxing the previous constraint, the intrinsic scatter also by varying 7', reaches 
its maximum for a DM homogeneous distribution (d=0) and its minimum for a 
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10 10.5 11 11.5 12 10 10.5 11 11.5 12 

log [M./MJ log [M./M e ] 




10 10.5 11 11.5 12 

log [M./MJ 



Figure 4: Intrinsic dispersion due to the logarithmic slope d of DM distribution in the three 
cases: d = (top-left), d = 0.5 (top-right), d = 0.897 (bottom-left) by 7' distribution frozen 
to that of the reference curve m° = 4. The mean dispersion values at the reference point 
Mb = Mg are given in Tab. 5. r becomes for 7^ = 2.103 when d=0.897 (see, text). See the 
on-line version for colours. 



2G 



10 10.5 11 11.5 12 10 10.5 11 11.5 12 

log [M./M B 1 log |M./M ] 

Figure 5: As in the previous Fig. 3 when the logarithmic slope d of DM density profile becomes 
d = (left) and d = 1 (right) (see, text). To be noted that the curve for m = m° = 4 remains 
the same whichever is the value of d; the others shift parallel by changing d if with the same 
m. See the on-line version for colours. 

cuspy halo with d — 1 (Fig. 5). That connection with DM cusp/core problem 
(see, Bindoni, 2008) is not new. In the CVMT indeed the tilt of galaxy Funda- 
mental Plane is obtained only with a cored DM distribution (d ~ 0.5) together 
with the explanation of many scaling relationships (see, SB9). On the contrary 
the observed tilt would disappear for cuspy halos with d = 1 . 

8. Kormendy relation 

Without ambition to treat here one of the main projections of FP, we under- 
line what is the (I e —r e ) relation we may immediately obtain from the mass-size 
relationship found. Indeed the CVMT leads to: 

M fl /Z~M£;a = ^ (50) 

From that, to get a Kormendy relation in the simplest way when: m = 4, 
d = 0.5, we use the corresponding linearized mass-size relation of the set (|46|) : 

r e ~ (L^f 755 ~ (7 e r e 2 ) - 944 (51) 
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which gives: 

h ~ r~ a94 ° (52) 

not too far (~ 24%) from what has been deduced from observations (see, 
D'Onofrio et al. 2006): 

I ~ r" 124 

1 e 1 e 

Many improvements may be done in this approach taking into account the 
two mass regimes and/or the curves at varying m instead of the linear approx- 
imation. 

To be noted that, even if the slope of mass-size relationship, i.e., the value 
p = 0.755, is the same for any d, nevertheless the result is depending on d via 
the tilt value of a. So, in the case of a cuspy halo, d = 1, i.e., corresponding 
to a not-tilted FP, the Kormendy's exponent of (|52|l turns to be —0.675 instead 
of —0.940 obtained with a cored halo of d = 0.5 which is far of about 45% from 
the observed value. 

9. The m-dispersion on other scaling laws 

The scaling laws for the three main quantities of FP: r e ,I e , a Q (j2~TI |2"51 |2T))) 
allow us to obtain not only the link between: Mb / L — a and L — a a (Faber- 
Jackson, FJ) but also to investigate the role of m-dispersion on these main 
scaling relationships and exspecially on the tilt of FP. Indeed, according to the 
CVMT the FJ-relation turns out to be: 

L ~ m(MV-i) aj y -'' (3 - d) (53) 

It is to be noted that both the exponents either of Mb or of m are directly 
related to the cosmological perturbation spectrum, via 7', and then they are 
depending on the amount of DM, but they are also function of how DM dis- 
tributes itself inside the halo, via the exponent d. So we have again to handle 
with a scaling law which, theoretically speaking, is not a straight line in the 
mass range of an observed sample of galaxies, but a curve having an intrinsic 
dispersion due to m. If we assume for a typical DM halo of Mo ~ 10 11 M Q , 
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that the typical mean value of 7' ~ 2 (ACDM scenario: n e ~ —2), then (|53p 
becomes: 

(d) a J(d) . (54) 



d = 0-=-l 



j = 0.44 ~ 



J = 2.7H-4 



If d = 0.5, we obtain: 

L ~ TO 016 a o 3 2 . 
It should be underlined that if d = 1 it turns out: 

without dependence on to. The range resulting for the a Q exponent J appears 
to be in good agreement with the value found by Faber (1987), 2.61 ± 0.08 (in 
B band) and with the significantly steeper slope found by Pahre et al. (1998), 
4.14 ± 0.22 (in K band). So the theory predicts an intrinsic different scatter in 
every FJ relationship due to the role of the factor m? which disappears as soon 



as the FP tilt disappears (d = 1 => a = 



l-d 
3-d 



0). 



Moreover the scaling laws ([2"7 ] US \ l2"9" ]l yield: 



M B /L - M% 



(3-d)- T ' 2 °t' 



which in turn gives for the gradient: 

dlog(M B /L) 2aj'(M D ) 



(55) 



(56) 



dlogo j'(Md) - 1 

It should be noted that, due to the proportionality of Mb with Md yn (in- 
cluding the fraction of DM inside at, see SB9), the same result holds for the 
log(Md yn /L) gradient. If again d — 0.5, then a = 0.2, and if 7' = 2, the gra- 
dient of eq.([56|) becomes 0.80, without regarding the factor m 04 . It is to be 
underlined that again for d — 1 together with the tilt the intrinsic dispersion 
in the relation (|56p disappears. A comparison test with the same gradient, as 
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derived from observations has been performed. Cappellari et al. (2006) found 
~ 0.8 for an observed sample of either fast rotators or non-rotating ETGs and 
SO galaxies. The fitted value by J0rgensen (1999) turns out to be 0.76 ± 0.08. 
Both are in good agreement with the theoretical result. 

On the contrary the ratio L/Mb which produces the tilt is totally indepen- 
dent of the cosmic perturbation spectrum and of the mass ratio m. It turns out 
to be dependent only on the dark matter density profile. Indeed, by definition 
it is: 

L/M B ~ I e r 2 jM B 

and by remembering the relationship between r e ~ a t ~ ao given by eq.([5]), the 
following functional dependence on the bright mass holds: 

L I Mb ~ M^^MjMs 1 ~ M B 1+t ^ (57) 

2 l'-(3-d) 

where the factor m ~<'( 3 - d ) introduced by I e is perfectly compensated by the 
factor hi V(3-d) which enters r\ (see (|27p and (J2SJ) whichever is the value of 
d. That is relevant into the explanation of the tightness of FP. 

Eq.(|57|) tells us how the ratio L/Mb loses its direct connection with the 
cosmology given by 7'. This confirmes the observation just made by Djorgovski 
(1992) that inside the FP as a whole there is a degeneracy in respect to the cos- 
mological process of galaxy formation which breaks into the projections of FP. 
Indeed the 7' springs out by considering what FP produces into the conjugate 
planes: i.e., the Faber- Jackson and the Kormendy relation. 

The above considerations are restricted to the effect of interaction energy 
between DM and star subsystem on the FP. Other physical ingredients induce 
surely a spread, such as the parameters of stellar populations, i.e. age and 
metallicity, rotational kinetic energy, radial orbital anisotropy and projection 
effects due to lack of spherical symmetry (e.g., Saglia et al., 1993). A detailed 
investigation on this point is outside the aim of the current paper. 
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10. Discussion 

Before conclusions we have to take briefly into account: i) the relationship of 
CVMT, in which two-component virial equilibrium holds conserving homology, 
while the weak homology appears in the FP literature (see, e.g., Bertin et al., 
2002); and ii) the consistency of CMTV interpretations in front of different new 
sets of data. 

i) It is important to distinguish between the dynamical effect the TCVM is 
able to offer in order to explain the tilt of FP without breaking the homology and 
the description of how the weak homology holds by adopting Sersic luminosity 
profiles. We mean: the dynamical theory tells us why the M/L ~ M a holds but 
it is not able to enter the details of how the light is distributed on the galaxy 
surface. On the other side, Sersic luminosity profiles are able to describe how is 
the best distribution of light (e.g., Kormendy et al., 2009) without explaining 
us why the total light L has to change by increasing M, i.e., L <~ M 1 ~ a . They 
appear as two complementary aspects of the main problem which is the change 
of light with mass. The description of Sersic does not enter the cause of the 
phenomenon, the explanation of the phenomenon does not enter the detailed 
description of its consequence. 

ii) We refer to the huge amount of database SDSS, published on the paper 
of Bernll. There are at least two very interesting outputs which spring out 
from the new relationships here derived by fitting the data. The first one is 
the correlation obtained for size and velocity dispersion a (see, Bernardi et al., 
2003a). From Fig. 2 of Bernll the correlation turns to be: r e ~ tr 15 in the 
mass range logM* = 10.5 = 11.3. Taking into account the Eqs. f([27f and (|29|) ) 
of the scaling laws, the prevision of the theory for the 0~ o exponent, in the same 
mass range, in the case of m = 4, disregarding the curvature of the mass-size 
relation, is: 2R/(1 — R) ~ 1.8. The second one is connected with the slope of 
the mass-size correlation once extended to the low mass regime. In the present 
paper we may evaluate the local slope down to (Mb/Mq) 9 - 60 . Disregarding the 
curvature it turns to be: p — 1/7' = 0.437 in the case m = 4 in agreement 
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with Bernll (Fig.l) where it has to be less than about 0.55. 
11. Conclusions 

By Clausius' virial maximum theory we have derived a — r e theoret- 
ical relationship for ETGs in a ACDM cosmological scenario which has been 
compared with that found by fitting the binned averaged data of PS96 (Tal09, 
Fig. 4). It is to be underlined that our theoretical slope changes with and 
increases with it. That is due to the increasing of the local slope in the cos- 
mological mass variance at increasing DM halo mass. The connection between 
cosmological scenario and luminous mass-size relationship is one of the best re- 
sult of the CVMT theory. Inside it many mass-size relationships are possible 
depending on the given Dark/Bright mass ratio m. The corresponding curves 
may be linearized in some fruitful numerical relations with a maximum scatter 
in respect to the right curves less than 1.5 times the rms observational errors 
(m = 1 -T- 4). Assuming the existence of two mass regimes, upper and lower 
10 11 M© according to Tal09, the agreement at lower regimes turns to be un- 
satisfactory, probably due to the crossing of different mass-size relationships at 
different m in the numerical fit of Tal09. At higher regimes the agreement turns 
to be fairly well. The difference between lower and higher regimes, underlined 
by Tal09, seems to be not so relevant from the theoretical point of view, due 
to the gradually increasing of the slope passing from one regime to the other 
without solution of continuity. Moreover, mass ratios m around 1—4 are more 
suitable to the galaxy sample used by Tal09, while values of about 10 appear to 
be ruled out, according to a ACDM scenario. 

The existence of individual galaxies with different values of m actually pro- 
duces a real scatter in the data we can see in the non-binned points (Fig. 4 of 
Tal09). This precious information is then locked inside the binning procedure. 
It is relevant that the amplitude of this dispersion is to connect with the dark 
matter inner distribution and more precisely that it increases or decreases going 
towards one cored or cuspy, respectively. 
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The last but not the least consideration has to be drawn about the capability 
of CVMT: the spread in the single galaxy data appears to be limited better 
from the lower part of the (M* — r e ) plane in respect to the higher one. The 
CVMT tells us that there has to exist a minimum value of m in the sample 
due to the presence of a sharper border on the side of lower r e . But what is 
its physical meaning? We know (Bindoni, 2008) that as soon as m decreases 
the ratio x = as I an increases towards 1. A limit mi exists and corresponds 
to as = old- After this limit the CV maximun can not exist because the 
equipartion of Clausius' virial energy is no longer possible. That limit value is 
mi < 1 and has to corresponds to the ZOE in the k — space rapresentation of 
FP. A reflex of ZOE is then also present inside the mass-size relationship with a 
'"Zone of Avoidance'" delimited by the mi — line (e.g., D'Onofrio et al., 2006). 

Many problems are still open in our approach. For example, the assumed 
homology requests to choose one common value for the King's concentration 
Cb and a unique CV maximum in the whole galaxy family to assure the same 
theoretical FP which implicitly underlies. But it is important to explore what 
happens by changing them. 

As general comment we can underline that the message which tells us this 
mass — size relationship is again that the true challenge is not to interpret in 
some way it but to understand why there are so many scaling relations of this 
kind for galaxies. In our opinion that is impossible to do whithout the cosmolog- 
ical framework because their existence is stricktly connected with how galaxies 
form. The Clausius' virial maximum theory included this link by nature. Based 
on classical physical grounds here revised the CVMT looks like very promising 
in order to understand the general features related to all astrophysical virialized 
structures and collected by Burstein et al. (1997) as cosmic metaplane. 
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